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^ ■ ABSTRACT 

^ \ We present a large-scale study of the X-ray properties and near-IR-to-radio spec- 

tral energy distributions (SEDs) of submillimetre galaxies (SMGs) detected at 1.1mm 
with the AzTEC instrument across a ^1.2 square degree area of the sky. Combin- 
ing deep 2 - 4 Ms Chandra data with Spitzer IRAC/MIPS and VLA data within 
> ■ the GOODS-N, GOODS-S and COSMOS fields, we find evidence for AGN activ- 

ity in '^14 percent of 271 AzTEC SMGs, '^28 percent considering only the two 
GOODS fields. Through X-ray spectral modelin g and multi-waveleng t h SE P fitting 
using Monte Carlo Markov Chain techniques to 'Sieben morgen et al.l ()2004f ) (AGN) 
and Efstathiou, Rowan-Robinson & Siebenmorgen (2000|) (starburst) templates, we 
. find that while star formation dominates the IR emission, with star formation rates 

' (SFRs) ~ 100 — 1000 Mq yr""'^, the X-ray emission for most sources is almost exclu- 

sively from obscured AGNs, with column densities in excess of 10^^ cm~^. Only for '^6 
percent of our sources do we find an X-ray-derived SFR consistent with NIR-to-radio 
SED derived SFRs. Inclusion of the X-ray luminosities as a prior to the NIR-to-radio 
SED effectively sets the AGN luminosity and SFR, preventing significant contribution 
from the AGN template. Our SED modeling further shows that the AGN and starburst 
I templates typically lack the required 1.1 mm emission necessary to match observations, 

" " arguing for an extended, cool dust component. The cross correlation function between 

the full samples of X-ray sources and SMGs in these fields does not indicate a strong 
correlation between the two populations at large scales, suggesting that SMGs and 
AGNs do not necessarily trace the same underlying large scale structure. Combined 
with the remaining X-ray-dim SMGs, this suggests that sub-mm bright sources may 
evolve along multiple tracks, with X-ray-detected SMGs representing transitionary 
objects between periods of high star formation and AGN activity while X-ray-faint 
SMGs represent a brief starburst phase of more normal galaxies. 

Key words: galaxies: active - galaxies: starburst - galaxies: high-redshift - submil- 
limetre: galaxies - X-ray: galaxies 
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1 INTRODUCTION 

Large blank-field surveys made at (sub-) millimetre wave- 
lengths have identified a large population o f bright, high- 



200C 



redsh i ft galaxies (e . 

20061 : iBertoldi et al] 1200 



Hughes et al.l 199^: ICoppin et al 



Perera et al] l2008l : IWei^ et al 



Scott et al. 2010l . and references therein). These 



sub- millimetre galaxies (SMGs) are cha racterized by high 
infrared (IR) lumino sities, >10^^ Lq (|Blain et al.l 12004 
IChapman et al.l |2005|). and a redshift distribution peak- 
ing around z ~2 l|Chapman et al.ll2005l V SMGs are there- 
fore believed to be the high-redshift analogs to local ultra- 
luminous IR galaxies (ULIRGs) and are possible progen- 
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itors of today's mass ive ellipticals fe.g. ISmail et al] |2004| : 
IChapman etaLll2005h . However, SMGs at z ~ 2 are more 
numerous than local ULIRGS by several orders of mag- 
nitude and lik ely dominate the tot a l IR luminosity den- 
sity a t z ~2 llLe Floc'h et al.1 l2005l : [ Perez-Gonzalez et alj 
I2OO5I : iHopkins et al.ll2010l ). The origin of these luminous, 
high-redshift sources is still under debate due, in part, 
to the low angular resolution at (sub-)millmetre wave- 
lengths of current instruments and the relative faintness 
of likely counterparts. Multi- wavelength and IR spectro- 
scopic follow-up s tudies of SMGs using Spitzer (see, for ex- 
ample, [M cncndc z-Delmestre et al] l2007l : iPope et al] l2008l : 
iNardini et all 12001 7 suggest that SMGs are largely dust- 
obscured starburst systems with star formation rates (SFRs) 
~1000 Mq yr~^. However, it is becoming increasingly ap- 
parent through the high X-ray detecti on rate of SMGs 
(~ 30-50%, see [Ale xander ct al. 2005a , b ^Laird et al.|[201ol : 
iGeorgantopoulos. Rovilos fc Comastrill201ll7 and SMG case 
studies (i.e., Tamura et alj 2010l ) that emission from active 
galactic nuclei (AGNs) may also be a crucial component to 
the energetic output of SMGs. 

The likely connection between starburst and AGN ac- 
tivity in SMGs is further supported by the concurrent 
nature of the cosmic SFR and black hole accretion wit h 
peaks at z ~2 (e.g.. iLe Floc'h et al.|[2005l : lMerlonill2004 ). 
Simulations of SMG formation in a merger-driven sce- 
nario also suggest that t he SMG phase preced es rapid 
growth of a central AGN (|Naravanan et al.1 12OI0I '). SMGs 
may therefore represent an important phase in galaxy evo- 
lution and may shed light on the origin of observed re- 
lations between AGN activity and stellar mass in local 
galaxies (i.e. the M-a relation; Ferrarese fc MerrittI I2OO0I : 
iGebhardt et~al]|2009l : [Cuhekin et aLll2009l ). One should be 
cautious, however, in extrapolating the starburst-AGN con- 
nection to the mo st extreme objects (i.e. radio-loud AGN, 
iDicken et al.ll201^ . and references therein) though such cases 
are a fundamentally different population of sources. Un- 
fortunately, while there are a multitude of methods for 
studying AGN and star formation, disentangling their rel- 
ative contributions to a galaxy's bolometric output re- 
mains challenging. Obtaining redshifts and other informa- 
tion via optical/ultra- violet imaging and spectroscopy is ex- 
ceptionally difficult a s SMGs are both distant and optically 
thick (see review by iBlain et alj|2002l ). IR spectroscopy of 
SMGs typically show strong polycyclic aromatic hydrocar- 
bon (PAH) features associated with star-forming regions, al- 
thoug h there are cases of power-law-like spectra indicative of 
AGN llMenendez-Delmestre et alj|2007l : ICoppin et alJbOOsI : 
iNardini et al.ll2008l : |Pope et al.l l2008^. Arguably, the best in- 
dicator for AGN activity is hard X-rays (>2 keV), which 
penetrate obscuring dust up to the Compton thick limit 
(neutral hydrogen column densities of Nh J^IO^* cm^^). 
X-ray detections are not uniquely attributable to AGN, 
however, as high SFRs may produce numerous high-mass 
X-ray binaries (HMXBs) that mimic the emission of low- 
luminosity AGNs. 

In the past decade there have been few studies that 
consider X-ray counterparts to SMGs for evidence of AGN 
activity, though t his num ber has expanded in recent years. 
[Alexander et al] l|2005al lbh (hereafter A05a,b) provide the 
earliest anal ysis by examining the Chandra counterparts 
to SCUBA (Holland et all \l99^ ) 850fim identified sources 



in the Great Observatories Origins Deep Survey (GOODS) 
North field. In their sample of 20 SMGs wit h radio and spec- 
troscopic redshift identifications taken from lChapman et al.l 
( 2005 ), they find that ~75 percent have X-ray properties 
consistent with obscured (Nh ^10'^'^ cm"'^) AGN activity. 
Accounting for SMGs without spectroscopic redshifts, they 
suggest that the true X-ray detection rate may be signif- 
icantly lower, >28 percent. However, the A05a,b sample 
may c ontain biases introduced through the 1 Chapman et al.l 
(j2005h SCUBA source catalog, which consists of obser- 
vations of known radio sources and low signal-to-noise 
(S/N<3.5cr) sources and thus may not be re presentative of 
the e ntire bright SMG population (see also, lYounger et alJ 
l200'it l. Furth er X-ray /SMG cou nterpart analysis has been 
provided by iLaird et al] (I2OIOI ') (hereafter LNPSIO), who 
find a ~45 percent X-ray detectio n rate to ra dio and/or 
Spttzer-identi&ed SCUBA sources l|Pope et alj [20061 with 
a ~20-29 percent AGN identification rate based on X-ray 
spectral modeling. LNPSIO find that the bolometric FIR 
emission is dominated by star formation in the majority of 
their sources (~85 percent) after including available Spitzer 
photometr y; consistent with A05a,b and other IR studies of 
SMGs (i.e. iMenendez-Delmestre et alj [20071 : IVahante et all 
I2OO7I : iMenendez-Delmestre et aLll2009l ). 

More recently, the stu dies of 

IGeorgantopoulos. Rovilos fc Comastril (120111) (h e reafter 
GRClll. iHiU fc Shanksl (|201lh and iBielbv et all (|2012l ) 
have utilized LABOCA d ata in the Exten ded Chandra 
Deep Field South (ECDFS, IWei/3 et al.|[2009h and William 
Herschel Deep Field. The analysis of GRCll is similar 
to that of LNPSIO who also find an AGN fraction of 
< 26 ± 9 percent with the mid-IR emission dominated by 
starburst activity, though the fraction of starburst-powered 
X-ray sources is lower t han es timated by LNPSI O . The 
works of Hill fc Shanki (1201 iD and Bi elbv et alj (|2012D 
consider a more statistical approach, utilizing the full 
catalogs rather than individual sources as in A05a,b, 
LNPSIO and GRCll, though find a similar SMG/X-ray 
detection rate (^^20 percent). They also find that obscured 
AGNs preferentially have greater sub-mm emission than 
unobscured AG Ns; a result confirm e d through ELVA 
observations by iHevwood erahl l|2012l) . iLutz et alj (I2OI0I ) 
find a similar relation in the ECDFS where the X-ray 
luminosity and absorbing column density for bright AGNs, 

is correlated with the 870/im 



L2 



> 10 



ergs s ^ 



fiux, implying a close connection to star formation. This 
assumes, however, that the X-ray emission is purely from 
the AGN while t he 870/jm flux is o nly from star formation. 
Furthermore, the lLutz et alj (I2OI0I ) study does not account 
for X-ray-bright SMGs, which may potentially bias the 
stacking results. 

To recap. X-ray studies to-date find that the AGN frac- 
tion of SMGs is in the range of ^^20-45 percent and that the 
bolometric IR luminosity of SMGs is dominated by star- 
bursts. 

In this work, we examine the identification rate and 
contribution of AGNs to the emission at various wavelength 
regimes in AzTEC SMGs. Our sample consists of Chandra 
X-ray counterparts to AzTEC 1.1 mm sources found in the 
GOODS-North, GOODS-South and COSMOS fields, pro- 
viding a total Chandra sky coverage of ~ 1.15 square degrees 
(~0.12, ~0.11 and ~0.92 square degrees, respectively) with 
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more than 2600 identified X-ray sources. This large sam- 
ple size will reduce any biases due to cosmic variance in 
previous studies. Furthermore, we do not base our sample 
selection and counterpart identification on prior source asso- 
ciation, thus removing any possible pre-identification bias. 
The available multi-wavelength photometry in these fields, 
including Spitzer IRAC and MIPS, will provide additional 
constraints on the AGN identification rate and contribution 
to the bolometric output of our sources. 

We begin with a description of the AzTEC and Chandra 
data and reduction procedures. We then detail our method 
for identifying X-ray counterparts to the AzTEC sources and 
subsequent multi-wavelength counterparts. Our analysis of 
the X-ray-identified AzTEC sources follows a two-pronged 
approach: (1) applying X-ray spectral models and SED tem- 
plates to the X-ray spectra and near-IR-to-radio SED, which 
will provide the basic information concerning the contribu- 
tion of AGN and star formation in each wavelength regime; 
and (2) linking the X-ray spectral fits to the near-IR-to- 
radio SED modeling, thus providing greater insight into 
the AGN/star formation connection. O ur SED fit ting dif- 
fers from typical SED analyses (e.g. iSerieant et al.i l2010) 
in that we employ a Monte Carlo Markov Chain (MCMC) 
technique. We close by comparing the implications of our 
work to those of previous X-ray/SMG results in addition to 
the X-ray/SMG cross-correlation relation. Additional analy- 
sis of our data, including source stacking and IR-optical-UV 
fitting, will be presented in future publications. 

Throughout this work, we assume a fiat ACDM cosmol- 
ogy with Ho = 70 km s~^ Mpc"\Ao = 0.73 and = 0.27. 



2 OBSERVATIONS AND DATA PROCESSING 
2.1 AzTEC: 1.1 mm Observations 



AzTEC (jWilson et al.ll2008l 'l is a 144-element bolometer ar- 
ray operating at 1.1mm an d installed on the 50m Large 
Millimetre Telescope (LMT: ISchloerb l'2008V Prior to its in- 
stallation on the LMT, AzTEC has performed several sci- 
ence observations on the James Clerk Maxwell Telescope 
(JCMT) and the Atacama Sub-millimetre Telescope Exper- 
iment (ASTE), including blank fields (namely GOODS-N, 
GOODS-S and COSMOS) and high redshift radio clusters. 
Here, we briefiy describe the AzTEC observations and 1.1 
mm source sample that will be used in our analysis. 

During the J CMT 2005-2006 observing campaign, 
iPerera et all (|2008l ) imaged a 21' x 15' area of the GOODS- 
N region. During the 2007 and 2008 obse rvation seasons on 
ASTE, AzTEC imaged both GOODS-S JScott et al.ll2010l) 
and t he one square degree area of COSMOS (lAretxaga et all 
I2OIII ). In reducing the raw time-steams for each set of obser- 
vations, an iterative technique using Principle Component 
Analysis (PGA) is used to filter out the at mospheric sig- 
nal th at dominates the raw observed data. iDownes et al.l 
(|201ll ) provides a discussion on correcting the PGA trans- 
fer function and lists revised catalogs for previously released 
AzTEC data. Here, we use the revised catalogs of Downes 
et al. for GOODS - N and GOODS-S; the COSMOS catalog 
of lAretxaga et al] (|201ll ) follows this prescription. The final 
AzTEC maps are constructed to have uniform coverage and 
sensitivity, providing a la rms of ~1.3 mjy in GOODS-N 



and COSMOS. The GOODS-S map reaches the confusion 
hmit of AzTEC on ASTE for a depth of (Icr) ~0.6 mJy 
rms. Sources are defined as peaks in the signal map with 
S/N^3.5(T, resulting in a total sample of 277 AzTEC sources 
(40, 48 and 189 in GOODS-N, GOODS-S and COSMOS, re- 
spectively) where <20 are expected to be false detections. 
Note, however, that the false detection rate is estimated for 
a S/N threshold of ~3.5cr and decreases rapidly for higher 
source S/N. For the following analysis, we use the full sam- 
ple of 277 AzTEC sources, applying no additional source- 
selection criteria. 

2.2 Chandra Observations 

The Chandra X-ray Observatory provides deep observa- 
tions of the GOODS-N, GOODS- S and COSMOS fields 
(for details on t he observations s ee [Alexander et al.l l2003l : 
ILuo et al.l I2OO8I : lElvis et al.l |2009| . respectively) with total 
exposure times of ~2Ms in each field. More recently, an ad- 
ditional ~2Ms has been added to GOODS-S with 31 addi- 
tional pointings; bringing the final integrated exposure time 
to ~4Ms (Xue et al., .20111 ). Due to the pointing strategy 
for COSMOS, effective exposur es only reach ~20 0ks for the 
inner ~0.5 sq. degree (see also lElvis et al.ll2009h . As a re- 
sult, the X-ray photon statistics in COSMOS are very poor, 
leading to weak constraints on the X-ray spectral properties 
(§ 3.1). This is somewhat offset by its larger area than the 
GOODS fields by allowing for more potential counterparts 
(§ 2.3). On the other hand, the deep 4Ms data in GOODS-S 
provides the greatest improvement to the counting statistics, 
and thus spectral modeling, to date; a valuable asset for po- 
tentially faint and highly obscured AGNs. All of the fields 
were imaged with the Advanced CCD Imaging Spectrome- 
ter Imaging array (ACIS-I), which is composed of four CCDs 
arranged in a 2x2 grid that operate together to provide a 
~ 17' X 17' field-of-view with sub-arcsecond resolution at 
the telescope aim-point, degrading with increasing off-axis 
distance. 

To ensure uniformity in our analysis, all observations 
were re-reduced using Chandra Interactive Analysis of Ob- 
servations (CIAO version 3.4) routines and custom routines 
developed for working with merg ed X-ray data s e ts; us - 
ing the published X-ray cata, l ogs of I Alexander et al.N2003l) . 
ILuo et al.l (|2008l ). IXue et all (|201ll ) and lElvis et al.l (|2009l ) 
would have required additional calibrations for compatibil- 
ity. Event files and exposure maps constructed in the 0.5-8.0 
keV energy range were made for all observations and then 
merged to produce final maps for the three fi elds. 

We use the source detection method of IWand l|2004h . 
with a false detection probability threshold of 10~^, to pro- 
duce X-ray source lists from the final images for cross- 
correlation with the AzTEC sample and spectral extraction. 
This detection method uses a wavelet analysis of the input 
images (in this case, the final merged X-ray images for each 
field) followed by a sliding-box map detection and maximum 
likelihood analysis for both source centroiding and optimal 
photometry. During the source detection, the X-ray maps 
are divided into different energy bands (i.e. 0.5-8.0 keV full 
band, 0.5-2.0 keV soft band and 2.0-8.0 keV hard band), 
resulting in a source catalog that includes all sources found 
in each energy band along with their respective count rates 
and positional uncertainties. The source detection process 
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also produces a list of source regions, which are defined as 
circular regions with radius equal to twice the 90 percent 
energy encircled fraction (defined according to the PSF at 
the source position). 

COSMOS poses a dilemma for source detection due to 
the blending of PSFs from the tiling of observations. To 
avoid this issue, we perform the X-ray source detection on 
the individual observations and then combine the resulting 
source lists into a final catalog. Derived parameters are re- 
calculated for each source using the final COSMOS map, 
with extraction radii determined from the smallest PSF 
corresponding to each source. Alternatively, one could sim- 
ply average the sub-catalogs to produce the final catalog; 
however, this may exclude X-ray counts present in an im- 
age where the source was not initially detected. Certainly, 
this method has difficulty in detecting the faintest sources 
present in COSMOS; nevertheless, this will not significantly 
influence our results given the already low depth of COS- 
MOS compared to the two GOODS fields. 

Combining the source lists from each field results in a 
total of 2630 X-ray sources available for our study. Individ- 
ually, there are 478, 526, and f626 sources in GOODS-N, 
GOODS-S, and COSMOS, respectively. Despite the differ- 
ences in data reduction and source detection, our source lists 
recover >90 percent of t hose from t he published cata l ogs o f 
[Alexander et al' ('2003VLu o et all (|2008l '), IXue et all (|201ll ') 
and lElvis et al. (2009). However, we miss many faint sources 
from the published catalogs due to our more stringent false 



detection threshold of 10 
catalogs. 



versus ~ 1-2x10 for the other 



2.3 Counterpart Candidates 

2.3.1 Chandra Counterparts 

The beam size of AzTEC on the JCMT and ASTE is 18" 
and 28" FWHM, respectively, making reliable X-ray coun- 
terpart identi f icatio n challenging. Following the method of 
IChapin et all (|20o3), we use a fixed search radius of 6" in 
GOODS-N and 10" in GOODS-S and COSMOS to find po- 
tential counterparts to the AzTEC sources. Our choice of 
10" in GOODS-S and COSMOS is consistent with a de- 
rived search radius fo r a source with S/N ^5.5 on the ASTE 
telescope according to llvison et al.l (|2007l ) and roughly corre- 
sponds to the average search radius for the AzTEC GOODS- 
S catalog JScott et al.l 2O10l). Simu lations in each field agree 
well with the Ivison et al.l |20o3) estimate and show that 
sources with S/N>3.5 are recovered within the respective 
search radii >85 percent of the time. Extending the search 
radius beyond our adopted value increases the number of 
X-ray counterparts; however, these additional X-ray sources 
are unlikely to be true counterparts (see below). 

As shown in Figure [TJ there is significant overlap be- 
tween the AzTEC and Chandra maps. Considering only the 
overlapping regions, our sample is limited to 271 (39, 47, and 
185 in GOODS-N, GOODS-S, and COSMOS, respectively) 
of the initial 277 AzTEC sources and 2229 (397, 429, 1403 
respectively) of the 2630 Chandra sources. Of the remain- 
ing 271 AzTEC sources, we find 38 with at least one X-ray 
counterpart (8, 16, and 14 for GOODS-N, GOODS-S, and 
COSMOS, respectively); 5 have 2 potential counterparts and 
1 has 3. For those sources with multiple potential Chandra 
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Figure 1. Chandra (solid black line) and AzTEC (dashed blue 
line) coverage regions for GOODS-N (upper), GOODS-S (middle) 
and COSMOS (lower). The AzTEC coverage given here corre- 
sponds to the 50 percent uniform coverage region used for source 
detection. Small circles with radii equal to the AzTEC beam- 
size (18" in GOODS-N and 28" in GOODS-S and COSMOS) are 
plotted at the AzTEC source positions. X-ray source positions 
are indicated hv the .small 'nhis' svmbols. 



X-ray Detections of SMGs 5 



counterparts, we treat each source individually and do not 
attempt to split the AzTEC flux as we have no prior in- 
formation on how it may be related to the potential X-ray 
sources. Overlapping spectral regions for these sources is not 
an issue as the uncertainty in the X-ray spectra is dominated 
by the low counting statistics. There are a total of 45 X-ray 
sources associated with the AzTEC sample, of which only 
2-3 are expected to be false identifications due to random 
alignments. Comparatively, the expected number of X-ray 
pairs for the entire sample of 271 AzTEC sources, assum- 
ing a purely random X-ray source population, is ~14. The 
AzTEC/X-ray identification rate is therefore ~14 percent, 
lower than estimates reported by A05a and LNPSIO due to 
the shallower X-ray depth of the COSMOS field; removing 
it increases the identification rate to ~28 percent. 

To assess the robustness of our X-ray counterpart iden- 
tifications, we compute the probability P of random associa- 
tion for a given AzTEC/X-ray pair given the search radii and 
X-ray source densities (2.97, 3.14 and 1.39 xlO~* arcsec"^ 
for GOODS-N, G OODS-S and COSMO S, respectively) us- 
ing the method of iDownes et al ] (|198^^. which corrects for 
the use of a finite search radius and fiux-limited source den- 
sity. The majority of the AzTEC/X-ray pairs (32/45) have 
P^0.05 which we define as a 'robust' counterpart, the re- 
maining AzTEC/X-ray pairs, with P=0.05-0.10, are 'tenta- 
tive' associations. Table [1] provides the list of the Chandra- 
detected AzTEC sources along with their relevant source 
properties and P values. 

Through this counterpart analysis, we are implicitly as- 
suming that the AzTEC and X-ray source populations are 
physically associated and that the two populations are not 
significantly clustered. If, on the other hand, the X-ray and 
SMC source populations are clustered, then we are more 
likely to falsely associate sources and misinterp ret the re- 
lation between AGN and starburst systems. ,Almaini et al.l 
(|2003l found evidence for correlation between Chandra and 
SCUBA 850^m source populations in the European Large 
Area ISO Survey (ELAIS) N2 field at the 4.3cr significance 
level and thus concluded that while they trace the same large 
scale structure, the AGN and starburst phases are not nec- 
essarily co-existent. Based on our cross-correlation analysis 
(see § 4.2), we find no evidence for significant correlation 
between deep Chandra and AzTEC source populations in 
general. 



(24 ^lm) SIMPLE0, GOOD^, and FIDElH data, in- 
cluding spectr oscopic /photometric redshift catalogs where 
available (e.g. *Barger ct al. 2003; B arger. Cowie fc Wand 



I2OO8I : ISantini et al. 2009; Silverma n et all 120101 '). Multi- 



2.3.2 Multi-wavelength Counterparts 

Thanks to the extensive multi-wavelength coverage in the 
GOODS and COSMOS fields, we are able to supplement 
the millimetre and X-ray data of our AzTEC sample 
with additional photometry and spectroscopic/photometric 
redshift^s from the GOODS and COSMOS public data 
sets. Accurate redshifts are the most crucial given the 
broad redshift distribution of SMGs and the sensitivity 
of X-ray spectral modeling to redshift (§ 3.1). Across the 
three fields, we u t ilize publicly available VLA (1.4 GHz; 
Miller et al.] I2OO8I : iKellermann et al.l I2OO8I : iMorrison et~al] 
2010l ). Svitzer IRAC (3.6, 4.5, 5.8, 8.0 ^m) and MIPS 



wavelength counterparts and redshifts for COSMOS were 
obtained by cro ss-referencing our detected sources with 
lElvis et alJ (|2009l ) and the COSMOS team's web-based data 
repositorylj In cross-referencing our AzTEC/X-ray sources 
with other catalogs, we use a search radius of 2", the av- 
erage X-ray positional uncertainty, centered on the X-ray 
counterparts. For each potential AzTEC/X-ray pair, we find 
no more than one potential counterpart in the VLA and 
Spitzer catalogs; these sources have been cross -checked with 
other AzTEC counter part publications (i.e. IChapin et al.l 
I2OO9I : lYun et ai] l2012t) and show excellent agreement. For 
reference, < 1 VLA./ Spitzer source is expected to be a mis- 
association due to random alignments over all three fields. 
For cases where we have IRAC but no MIPS identifications, 
we estimate a 5(t MIPS flux upper limit through the photo- 
metric error of the MIPS source nearest to the IRAC posi- 
tion. A complete catalog of the multi-wavelength photome- 
try and redshift data for our sample is given in Table [2] 



3 ANALYSIS 

With our sample of X-ray selected AzTEC sources in hand, 
we now examine their physical properties through a variety 
of methods. We start with modeling of the X-ray spectra. 



3.1 X-ray Spectral Modeling 

X-ray sources with L2.0-10.0kcV 10"*^ ergs s~^ are generally 
believed to be powered almost exclusively by AGN with ab- 
sorption due to modest amounts of dust and gas within the 
host galaxy. A05b showed that X-ray-identified SMGs are 
predominately heavily obscured, possibly even to the Comp- 
ton thick limit with column densities of Nh ^ 10^^ cm~^. 
For the most extreme cases of obscuration, a buried AGN 
may only be visible in light scattered off of the obscuring 
torus. Alternatively, if SMGs are powered by a high rate 
of star formation, then the observed X-ray emission could 
result from the stellar population, powered by numerous 
high- mass X-ray binaries (HMXB) . For comparison, a typi- 
cal SMG with SFR in the range of 100-1000 M© yr"^ would 
produce an X-ray source with 2.0-10 .0 keV luminosity of 
~ 1041-42 ergs s"^ (|Persic et al.ll2004 hereafter P04). 

For our sample of AzTEC/X-ray sources, we first ex- 
tract their source and local background spectrum in the 0.5- 
8.0 keV observed energy range using the region files defined 
from our source detection (see § 2.2). Note that background 
spectra are taken from source-removed event files to avoid 
contamination from nearby s ources. The spectra are fitted 
in the xspec (version 12.4.0, !Arnaud"l996l, |2003| ) software 
package using the C-statistic (|Cash 1979: due to the low 
photon counts in many of the spectra (see Table [l| . In or- 
der to improve the counting statistics within each bin, we 



1 http:/ /www. astro. yale.edu/dokkum/simpIe/ 

2 http:/ /www. stsci.edu/science/goods/ 

^ http:/ /ssc. spitzer. caltech.edu/legacy/abs/dickinson2. html 
* http://irsa.ipac.caltech.edu/Missions/cosmos.html 
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Table 1. Chandra identifications of AzTEC sources in GOODS-N, GOODS-S and COSMOS. Errors are given at the Icr confidence level. 
Col.(l): AzTEC source ID prefixed by field (i.e. AzGN24 for source 24 in the AzTEC GOODS-N catalog). Col.(2): Chandra ID following 
lAU standards. Col. (3): Positional offset between AzTEC and Chandra sources. Errors are derived from Chandra positional uncertainty. 
Col. (4): Chandra 0.5-8.0 keV full band count rate. Cols. (5): Total counts within the source regions as defined from our X-ray source 
detection. Col. (6): Estim ated background counts within th e source regions. Col. (7): deboosted AzTEC source flux (see section 3.5 of 
lAustermann et al.|[201ol and section 6.2 of lScott et "al]|2010l) . Col.(8): Probability P of the Chandra source being a random association. 
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t Source also detected in LNPSIO. 



have re-binned the spectra to fixed width spectral channels 
of ~43.8 eV. 

In fitting the X-ray spectra, we consider two differ- 
ent classes of spectral models: (1) an intrinsically absorbed 
power- law, indicative of AGN; and (2) a stellar model based 
on HMXB emission including intrinsic absorption. These 
models are designed to be simple, yet physically meaningful, 
representations of the X-ray emission. For comparison with 



previous works, we also consider a simple power-law with 
only Galactic absorption, represented by the xspec model 
pha(po), to measure the effective photon index Tsff- As 
the C-statisti c itself is no t a measure of the "goodness-of-fit" 
(see, however. ILucv|2 00(?1 . we use the xspec GOODNESS com- 
mand for comparing the different spectral models (§ 3.1.3). 
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Table 2. VLA, Spitzer IRAC/MIPS and redshift information for the X-ray identified AzTEC sources. Spectroscopic and piiotometric 
redsliift information for the AzTEC/X-ray sources was taken, primarily, from publicly available redshift catalogs (see § 2.3.2 for details). 
MIPS upper limits are estimated from the 5cr upper limit of a detected MIPS source nearest the AzTEC/X-ray position (§ 2.3.2). Errors 
are given at the la confidence level. 
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1, 


.033 
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1, 


.101 


1.330 
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58.2+0.3 
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61.0+0.5 


47.2+1.1 
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0.910 
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J100033, 


,75+014906.3 




















AzC17 


J100055, 


,34+023441.1 


78+12 


1390+20 


99.7+0.2 


166.1+0.4 


254.9+1.2 


407.6+3.0 


1, 


.404 


1.410 


AzC147 
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,46+015718.1 




80+10 


36.6+0.2 


35.7+0.3 


29.3+1.1 


29.3+2.3 






1.230 


AzClOS 


J100116, 


,15+023606.9 




520+60 


128.1+0.2 


140.6+0.4 


162.7+1.1 


188.7+2.5 


0, 


.959 


0.950 


AzC85 


J100139, 


,73+022548.5 


549+12 


180+20 


1100.3+2.3 


780.2+2.0 


510.3+2.1 


346.1+3.0 


0, 


.124 


0.120 


AzGll 


J100141, 


,02+020404.8 




210+20 


11.2+0.1 


16.3+0.2 


26.5+1.0 


40.9+2.3 









t The photometric redshift was adopte d for J033215. 32-27 5037. 6 following cross-catalog comparison with GOODS-MUSIC 

ijSantini et ahllioogl ) and additional analysis. 



3.1.1 Model A: Absorbed Power-Law 

Our first model provides a simple parametrization of tlie X- 
ray emission from an AGN, represented by a single power- 
law. The model includes the effects of both (Milky Way) 
Galactic and intrinsic absorption and is represented by the 
XSPEC model pha(zpha(po)). The X-ray spectra is thus de- 
fined by the intrinsic absorption, Nh, and photon index, F. 
As these values can be strongly correlated for weak sources, 
we chose to fix the photon index to F = 1.8, typical for 



unob scured AGNs (i.e. iNandra fc Poundsl[l993 : iTozzi et all 
I2OO6I ). The model (hereafter Model A) thus represents a typi- 
cal AGN and provides an estimate of the level of obscuration 
present in our X-ray-identified SMGs. 



3.1.2 Model B: A bsorbed HMXB 

Our second model (Model B) is developed for emission due 
to star formation and is based on the HMXB X-ray spectral 
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Model A 




Rest-frame Energy (keV) 



Figure 2. Comparison of the X-ray spectral Models A (solid) and 
B (dot-dasiied) normalized at ~10 keV. The models arc shown 
for fiducial column densities of lO^'^ cm~^ (top) and 10^^ cm~^ 
(bottom). The shaded region indicates the effective rest-frame 
energies sampled by the 0.5-8.0 keV observed spectrum of a source 
at z ~ 2. 

model of iPersic fc Rephaelil (|2002l ). In summary, the X-ray 
emission from HMXBs can be expressed as a broken power- 
law of the form 

where r aee=l-2 (typical of bright, ac cretion powered X-ray 
sources; I White. Swank fc Holtl [l98^ 'l with a cutoff energy 
of Ec ~20 keV and e-folding energy of ep ^^12 keV. Ide- 
ally, when constructing a spectral model for stellar pro- 
cesses, we should also include contributions from low mass 
X-ray binaries (LMXBs) and supernovae. However, super- 
novae contribute little to > 2keV rest-frame flux compared 
to HMXBs. While LMXBs may contribute a considerable 
fraction of the hard X-ray flux, the low mass stellar com- 
panion typically has not had time to evolve off the main 
sequence and fill its Roche lobe by z ~l-2. For sources in 
our sampl e with z < 1, we may st ill use the HMXB-SFR 
relation as iPersic fc Rephaelil l|2007l ) showed that for mod- 
erate to high SFRs (SFRs>50 Mq yr"^) the X-ray-SFR rela- 
tion is similar to the HMXB-SFR relation. Our stellar spec- 



tral model therefore consists of only the HMXB emission, 
which is absorbed by both Galactic and intrinsic material 
(pha(zpha(hmxb)) in xspec, where the model hmxb is de- 
fined as given above). We include intrinsic obscuration in 
Model B, since it is clear from multi-wavelength evidence 
that SMGs are heavily dust-enshrouded systems. With Vacc 
fixed to 1.2, we are left with only the intrinsic obscuration 
and normalization to vary between spectra, similar to Model 
A. 

As shown in Figure [2l there are immediate differences 
in the spectral shapes of our adopted models. Both models 
appear similar at low energies; however, the difference in 
spectral slopes, as well as the exponential cut-off in Model 
B, are apparent for higher energies. For high obscuration and 
low count spectra, it is difficult to distinguish between Model 
A and B (§ 3.1.3). However, the derived Nh values will vary 
according to the power-law spectral slope. Additionally, we 
can compare the X-ray-derived SFRs of Model B with those 
obtained through our NIR-to-radio SED modeling (§ 3.2). 

3.1.3 Application of X-ray Spectral Models 

We now apply our set of spectral models to the X-ray identi- 
fied AzTEC SMGs. To correctly fit the intrinsic absorption, 
which has a strong energy dependence through the photo- 
electric cross-section, we require accurate source redshift in- 
formation. This limits us to 32 out of our original sample 
of 45 X-ray sources (~63 percent), including 5 sources in 
GOODS-N, 14 in GOODS-S, and 13 in COSMOS. We favor 
the spectroscopic redshift, whenever available, over the pho- 
tometric redshift. Milky Way absorption values of 1.5, 0.9, 
and 2.5x10^" cm~^ are included for the spectra, depending 
on whether they were taken in GOODS-N, GOODS-S, or 
COSMOS, respectively. The best-fit parameters for each set 
of models, as well as their C-statistic values and associated 
rest-frame, absorption corrected 2.0-10.0 keV luminosities, 
are given in Table O As a simple check, we have compared 
our derived l uminosities with those of previously pu blished 
catalogs (i.e. [Alexander et al. l l2003l : lTozzi et al.ll2006l ) which 
correlate well with our results. 

In order to determine which of our sets of models offer 
the best fit to the X-ray spectra, we run 2000 Monte Carlo 
simulations through the GOODNESS command in xspec, 
which provides the percentage of simulations that have a 
C-statistic lower than the observed spectrum. The best-fit 
spectral models, the ones providing the lowest goodness frac- 
tion, have been highlighted in boldface in Table [S] As one 
might expect, the models with the lowest C-statistics tend to 
also provide the lowest goodness fractions, indicating a very 
high probability that the observed spectrum can be charac- 
terized by the best-fit model. Models A and B often show 
very similar C-statistics, which leads to only a few percent 
difference in their goodness fractions. These differences are 
not statistically significant based on 10000 fakeit simulated 
fits using an intrinsically absorbed r=1.8 power-law as the 
template spectrum. 

We find that ~53 percent (17/32) of the AzTEC/X- 
ray sources have X-ray spectra that immediately favors 
an AGN origin. Of these, ~70 percent show evidence for 
heavy obscuration with Nh > 10^^ cm~^. Regardless of 
the best-fit spectral model, the majority of AzTEC/X-ray 
sources (22/32) have 2.0 to 10.0 keV rest-frame luminosi- 
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Table 3. X-ray spectral fits to identified AzTEC/X-ray sources. Spectral models used are: Galactic dust- and intrinsically-absorbed 
AGN power-law (pha(zpha(po)), Model A) and Galactic dust- and intrinsically-absorbed power-law with an exponential cut-off relating 
to emission from HMXBs (pha(zpha(hmxb)), Model B). Models that offer the best fit to the X-ray spectra based on our simulations are 
emphasized in bold. The relevant parameters given here are the intrinsic neutral hydrogen column density (Njj in 10^^ cm~^); absorption 
corrected, rest-frame X-ray luminosity in the 2.0-10.0 keV energy band (Lx in 10^'^ ergs s~^) and X-ray derived SFR (SFRx in 1000 
M0 yr-l) for Model B assuming the P04 relation. Errors are given at the 90 percent confidence level. 



Chandra ID Ve f f Model A Model B 









Nh 


Lx 


C-stat 


Nh 


Lx 


SFRx 


C-stat 


J123616 


08-1-621514.1 


0.98+n?a 
— 0.28 


16.46+^-^;^ 
—6.20 


4.40 


159.3 


6. 70+"? 

— 4.05 


2.60 


26.0 


161.0 


J123635 


86-1-620707.8 


-0 56+°-^? 


97 94+25 !? 


8.89 


212.7 


74 19+25--'-^ 


4.51 


45.1 


213.3 


J123711 


32-1-621331.1 




9 72+2"'".;?^ 
^- ' -6.95 


0.92 


189.0 


< 16.28 


0.62 


6.2 


186.7 


J123711 


98-1-621325.8 


-0 41+2 f 


01 ■'JU_24.30 

2 10+°-25 
-^•■'■"-0.23 


2.70 


199.2 


38 31+?I-°S 


1.37 


13.7 


198.8 


T1 9*^71 

d l^O 1 -LU 




-, ,--1-0.05 
^■^'-0.06 


9.95 


236.6 




8.76 


87.6 


239.4 


J033158 


25-274458.8 


U.»»_Q gg 


< 1.93 


0.07 


177.3 


< 0.83 


0.08 


0.8 


173.8 


J033204 


48-274643.3 


1 QQ-l-0.24 


i.oo+,V9? 

■^■■""-0.82 


1.08 


187.0 


< 0.37 


1.07 


10.7 


186.1 


1033207 


12-275128.6 


74+;: 2-7^ 

^- '^-0.75 


< 14.0 


0.04 


189.1 


< 12.3 


0.04 


0.4 


189.2 


J033209 


71-274249.0 


2 22+°-^9 
^•^^-0.28 


< 0.06 


0.25 


199.4 


< 0.04 


0.34 


3.4 


238.4 


J033212 


23-274620.9 


85+°-^^ 


'-'■■-'^-0.68 


1.00 


189.8 


1 53+°-^8 

J..OO_Q go 


0.88 


8.8 


186.6 


1033215 


32-275037 fi 




"■°-'-0.42 


0.01 


204.2 


Q4+0-35 


0.01 


0.1 


205.7 


J033222 


1 7-27481 1 6 


u.oo_o 28 


39 19+14.60 
o».ia_jo 4Q 


3.11 


181.8 


23 67+1° **^ 
■'■J-O' _g.72 


1.66 


16.6 


182.1 


J033222 


56-27481 5 


-0 43+0-49 

^•^'-'-0.42 


Q4 1 1 +55.73 


3.51 


182.9 


55 83+45-48 


1.49 


14.9 


180.6 


1033234 


78-275534 


1 06+°-^* 


43+0-27 

"■^■-'-0.19 


4.0e-4 


167.9 


< 0.45 


5.0e-4 


5.0e-3 


167.7 


J033235 


18-275215.7 


64+0-39 


5 17+4.38 

-2.08 


0.12 


179.7 


q 1 7+3.36 
•^--■-'-1.84 


0.10 


1.0 


180.9 


J033238 


01-274401.2 


1 7--I-1.00 
'-0.88 


< 5.53 


0.14 


235.8 


< 3.95 


0.14 


1.4 


237.6 


J033244 


02-274635.9 


^•"^-0.20 


< 0.96 


3.64 


181.6 


< 0.26 


3.15 


31.5 


230.3 


J033246 


83-275120.9 


95+0-52 


4.32t«^ 


0.20 


209.2 


< 5.64 


0.17 


1.7 


208.9 


J033302 


94-275146.9 






14.37 


175.5 


< 7.89 


8.69 


86.9 


182.4 


J095905 


05-1-022156.4 


0.98ii;|° 


< 78.54 


4.92 


47.1 


< 60.47 


2.72 


27.2 


47.6 


J095929 


70-1-021706.4 




< 4.14 


0.83 


92.7 


< 3.28 


0.93 


9.3 


91.8 


J095945 


15-f023021.1 


n t;7+0.58 
^•^'-0.59 


< 1.01 


0.41 


140.0 


< 0.95 


0.57 


5.7 


139.7 


J095953 


85-1-021853.6 


5.56+^?« 


0.79 


103.2 


q 07 + 3.66 

■^•-^'-2.56 


0.66 


6.6 


103.7 


J095959 


96-1-020633.1 


52+0-66 
^•3^-0.74 


5.5314;?? 


0.39 


79.0 


q 71+3.92 

'^-2.65 


0.33 


3.3 


80.0 


J100003 


73-1-020206.4 


i.ool?"^ 


< 124.89 


6.16 


141.2 


< 82.08 


2.52 


25.2 


140.8 


J100006 


11-1-015239.2 


-I 7--I-0.70 


< 1.48 


2.25 


115.6 


< 0.83 


2.26 


22.6 


118.3 


J100006 


55-1-023259.3 


1 26+°-'^* 


< 2.98 


0.82 


87.1 


< 1.56 


0.79 


7.9 


86.4 


J100033 


61-1-014902.0 


, c-y-l-O.SO 
-'-•3'-0.45 


< 0.68 


0.75 


104.8 


< 0.29 


0.95 


9.5 


107.5 


J100055 


34-1-023441.1 


1 ccr-l-0.19 
l-SO-0.19 


< 0.44 


24.88 


159.2 


< 0.10 


26.84 


268.4 


181.6 


J100107 


46-1-015718.1 


1 5Q+0.79 

j^-oy-0.60 




1.46 


148.5 


< 1.50 


1.49 


14.9 


149.9 


J100116 


15-1-023606.9 


1 79-1-0.60 
^•'^-0.56 


0.20lJ;?i 


5.84 


100.3 


< 0.69 


6.72 


67.1 


102.6 


J100139 


73-1-022548.5 


q 90-1-0.79 
•^•^'-'-0.71 


< 0.05 


0.01 


81.3 


< 0.04 


0.02 


0.2 


97.5 



ties of >1043 ergs s~^, heavily favoring an AGN interpreta- 
tion. Note that the derived luminosities are sensitive to the 
choice of the X-ray model. For those AzTEC/X-ray sources 
that favor the starburst model Model B, we use the X-ray 
luminosity to SFR relation of P04 to estimate a SFR, as- 
suming no contribution from a buried AGN. There is some 
uncertainty in the exac t form of the X-ray-to-SFR scaUng re - 
lation as discussed by iMineo. Gilfanov. fc SunvaevI (|201l[ ) ; 
however, many of these relations consider local, low SFR 
(< 10 Mq yr~^) sources during their construction. As we 
are concerned with potentially high SFRs, we favor the P04 
and iPersic fc Rephach (2 007) SFR-X-ray scaling relations; 
using the Ranalli et al.l ( 20031 ') relation, or similar, would 
decrease the estimated SFRs by a factor of ~2-5. The high 
X-ray luminosities would require very strong SFRs on the 
order of > 10'' — 10* M© yr~^, which is pushing the lim- 
its for typical SMGs. However, there are 5 sources with 
Lx < 10*^ ergs s~^ which are candidates to be starburst 
powered X-ray sources. These sources account for ~16 per- 



cent of our X-ray-identified SMG sample; consistent with the 
starburst-powered fraction of LNPSIO (~17±6 percent). We 
caution, however, that this does not necessarily imply that 
their X-ray emission is dominated by star formation (see 
§ 3.2.2, Tables |4]). These results thus show that the bulk of 
the X-ray emission from our SMG sample is predominately 
produced by obscured AGNs. 



3.2 NIR-to-Radio SED Modeling 

For an alternative view of the AGN and star formation con- 
tributions, we now examine the near-IR-to-radio SEDs of the 
AzTEC/X-ray sources. To be luminous at (sub-)millimetre 
wavelengths, a source must contain dust heated to T~30K 
(Chaoman ct al. 2005; Podc ct al. 2006) through some cen- 
tral engine. While it is possible to have (sub-)mm emis- 
sion due to synchrotron processes from radio-loud AGN 
(e.g., IVieira et al.l |2010[ ). the corresponding radio fluxes 
would have to be significantly larger (on order 1-100 mjy; 
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Ide Zotti et"ai1l2010l : IVieira et al1l2010l 'l than those observed 
for our AzTEC/X-ray sources, which range from 0.02 to 0.65 
mjy (Table [2)|. The required dust heating must then be ac- 
compUshed either by star formation, AGN activity, or some 
combination of the two. 

For our SEP modeUng, we consider the t em plate s 



of 
and 



ipiates 

fstathiou. Rowan-Robinson fc SiebenmorgenI (|200GI ) 
Siebenmorgen et al.l ( 20041 ') to parametrize emission 



from a starburst (SB) and AGN component, respectively. 
This selection of templates is widely used in the liter- 
ature and has shown to provide reasonable results to 
similar class es of sources over the NIR-to-mm wavelength 
regime fi.e. lEfstathiou. Rowan-Robinson fc SiebenmorgerJ 
2000l: ISiebenmorgen et all l2004l: iMeng et all I2OI0I: 



Serraetal.1 I2OIII : lYounger fc Hopkins! I2OIII : lYun et all 
201I and references therein). For this work, we favor the 



Siebenmorgen et al.l (|2004h AGN models as opposed to 



torus models as we are more interested in the integrated 
AGN host properties rather than the centralized nuclear 
region. Additionally, these models are built from basic 
radiative transfer models, incorporating relevant dust 
emission/absorption physics, with simple parametrizations 
comparable to the SB models. 

In order to estimate the total SED, we apply a simple 
linear combination of the two template sets. Since this ap- 
proach may introduce strong template-parameter degenera- 
cies into our summed SEDs, we use a Monte Carlo Markov 
Chain (MCMC) technique to perform the fitting. While 
computationally slower compared to direct maximum likeli- 
hood (least squares) fitting, MCMC has the advantage that 
the full set of posterior parameter distributions are returned 
- allowing for direct inspection of the posteriors for degen- 
eracies that may bias our interpretations of the fits (see Fig- 
ure |4]). The full details of this method will be presented 
in Johnson et al. (in prep). Here, we briefiy describe the 
adopted models and their implications on the AzTEC/X- 
ray source population. 



3.2.1 SED Models and Fitting 

Before applying the SED templates to the observed SEDs, 
it is helpful to have an understanding of how the templates 
parametrize the underlying physics and resulting IR emis- 
sion. In the lEfstathiou. Rowan-Robinson fc SiebenmorgenI 
l|2000 ^ templates, emission from a dusty starburst is traced 
from a single star forming CMC with the cloud optical depth 
iji,) and starburst age setting the overall shape of the SED. 
Specifically, Ti, controls the strength of the PAH and sili- 
cate features, while older starburst ages shift the IR peak to 
longer wavelengths. A normalization factor is then required 
to scale the emission from a single CMC to the full system. 
This normali zation is comparable to the SFR at the onset of 
the bu rst as lEfstathiou. Rowan-Robinson fc SiebenmorgenI 
l|2000l ) assume an exponentially decaying SFR history of the 
form 



SFRit) « SFR{0)e 



-t/20 Myr 



(2) 



where t is the SB age; SFR estimates obtained this way 
are approximately 2-3 times lower than more trad i tiona l 
FIR SFR indicators - for example, the iKennicuttI (|l998l ') 
relation. The AGN models are described by a single cen- 



tral illuminating source with intrinsic luminosity L sur- 
rounded by a spherical dust distribution of size R and 
the visual extinction (Ay). The dust distribution, temper- 
ature, strength of absorption/emission lines, etc. are ad- 
justed through a combination o f the size and visual e xtinc - 
tion. It should be noted that the lSiebenmorgen et al.l (|2004f ) 
AGN templates make a number of simplifications compared 
to alternative AGN models. Modern AGN templat es (e.g. 



Fritz. Franceshini fc Hatziminiaog lou 2006 ; Nenkova et al.l 
20081 . and references therein) consider the AGN to be sur- 



rounded by torus, generally composed of a clumpy mate- 
rial, whose geometry flares outward. This geometry nat- 
urally falls in line with the standard AGN unifled model 
where looking through the torus results in Type 2 (ob- 
scured) AGN while Type 1 (uno bscured) AGN are produced 
from 'face-on' observations. The lSiebenmorgen et al.l (|2004l ) 
models obviously lack the asymmetry and clumpy distri- 
bution of the traditional AGN torus but are able to recre- 
ate the same effects; Siebenmorgen et al. comments that it 
is the dust mass and distance from the source (set by Av 
and R) that are most important. Though torus geometries 
may extend to the kpc scale (e.g. .Gr anato & P anose 199 ^; 
Frit z. Franceshini fc HatziminiaoglorJ bood : iNenkova et al.l 
1200^ and can produce significant cold dust emission, they 
lack the dust intrinsic to the host galaxy, whose geometry 
extends well beyond that of a nuclear torus. Given that the 
photometr y for our sample can not re solve our sources, we 
believe the ISiebenmorgen et al.l (|2004| ) models to be better 
representative of galactic emission resulting from an AGN 
than the traditional torus models. 

Using the above models, we have three to six free pa- 
rameters with 6-7 available SED data points per AzTEC/X- 
ray source. In our fitting, we are able to predict an X-ray 
lu minosity from the FIR luminosity/SFR using the relations 
of lMarconi et al.1 (120041 ) and P04 for the AGN and SB mod- 
els, respectively. This allows us to then use the observed X- 
ray luminosities in Table |3] as an additional prior to the fits. 
As the SB models do not account fo r any radio emission, we 
emplo y the radio-FIR correlation of lYun. Reddv fc Condoiil 
(|200lh to add a radio 'tail' to the templates. Note, how- 
ever, that this may still pose some uncertainty when com- 
bining templates as there is scatter in this relationship (e.g. 
ICarilli fc Yuiill2000l : IChapman eral]|2005h and it does not 
predict any radio emission resulting from an AGN. 

In fitting the near-IR-to-radio SEDs of our X-ray- 
detected SMGs, we consider two combinations of the SED 
templates: (1) AGN and SB templates including the ob- 
served X-ray luminosity and X-ray-absorbing column den- 
sity as priors to the AGN luminosity, SB SFR and AGN Av 
and (2) SB only without the additional X-ray constraints. 
The first set of models serves to estimate the AGN contri- 
bution to the bolometric and 1.1 mm emission. The SB only 
fits provide a measure of the necessity of the AGN tem- 
plates. The X-ray luminosity prior had to be excluded for 
these fits as their inclusion produced unreasonable results 
(see § 3.2.2). Tables U and |B] and Figure [3] show the results 
of our MCMC fitting technique to our AzTEC/X-ray sam- 
ple. For each set of best-fit parameters, we calculate the log 
of the likelihood, ln(L); higher values of ln(L) indicate a 
higher probability that the data is consistent with the best- 
fit model. 
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Table 4. Best-fit parameters for the composite AGN+SB models based on the broadband photometry of AzTEC/X-ray sources. The 
predicted X-ray luminosities are compared to those derived from the X-ray spectral modeling (§ 3.1.4) to provide additional weights in 
calculating the likelihoods. Errors are given at the Icr confidence level after marginalizing over all other free parameters in the fitted 
templates. Col.{l): Chandra Source ID. Col. (2), (3), and (4): AGN template galaxy outer radius, intrinsic luminosity and visual extinction. 
Col. (5), (6), and (7): SB template normalization, age, and optical depth. 
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J123616.08-I-621514.1 
J123635.86-I-620707.8 
J123711.32-I-621331.1 
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J095905.05-I-022156.4 
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J100006.11-I-015239.2 
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J100055.34-I-023441.1 
J100107.46-I-015718.1 
J100116.15-I-023606.9 
J100139.73-I-022548.5 
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3.2.2 SED Fitting Results 

Figure [3] shows that our method is able to produce 
reasonable fits to the AzTEC/X-ray sources; the source 
J033234. 78-274815.0 was excluded as it has no discernible 
IRAC/MIPS counterpart despite having a spectroscopic red- 
shift (see Table [2]). While the majority (~87 percent) of our 
sources can be fit using the SB templates alone, they typi- 
cally under-predict the 1.1 mm emission, recovering on aver- 
age ~30-38 percent of the observed flux. Including the AGN 
models helps to slightly increase the model fluxes and are 
generally required to match the X-ray luminosity prior but 
are still unable to match the mm- wavelength observations; 
contributing little, if at all, to the bolometric luminosities 
and observed 1.1 mm flux. In some cases (e.g. J033212.23- 
274620.9), the AGN and SB templates appear very simi- 
lar in the flnal fit. This likely results from the similarities 
in the dust treatm ent and radiative transfer i n the tem- 
plates as noted by ISiebenmorgen fc Efstathioul (|2005h for 



effectively identical template parameters (i.e. dust content, 
optical depth and intrinsic luminosity; see also their figure 
4). In many cases, the fit values for AGN R and Ay are 
rather poor and show a large range in acceptable values. This 
effect stems from our use of the X-ray luminosity prior which 
effectively sets the AGN bolometric luminosity, preventing 
any additional AGN contribution to the bolometric SED and 
thus leading to unconstrained R and Ay (see also Figure [5]). 
The fact that Ay had such poor constraints prompted us to 
include the X-ray column density to avoid over-estimating 
the dust content. By virtue of our MCMC technique, we may 
readily identify any degeneracies between the AGN and SB 
template sets; however, Figure |3] shows that there are no 
large parameter-parameter degeneracies, although some pa- 
rameters are not very well constrained. This is particularly 
the case for AGN R and Ay as mentioned above. 

Due to the under-prediction of the 1.1 mm flux, we are 
biased to under estimate the total IR luminosity such that 
the values reported in Table |4] are more likely to be lower 
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Table 5. Continuation of Table l4l containing the derived physical properties based on the fitted parameters. Col.(l): Chandra Source 
ID. Col. (2): Total model-derived, rest-frame bolometric IR luminosity from the AGN and SB templates over the wavelength range 
~0.001-1500/.tm. Col. (3): SFR derived from Eqn. 2 (§ 3.2.1) and Cols. 5 & 6 of TableH Col. (4); SED derived X-ray luminosity Col.(5): 
Fractional contribution of AGN template to total model emission at 1.1mm. Col. (6): Fractional contribution of model to observed 1.1mm 
flux. Col. (7): ln(L) of best-flt parameters. 





IR 




X-ray 








Chandra ID 


Lum. 

IO^^Lq 


SFR 

Mq yr-1 


Lum. 

lO^^ergs s-i 


fAGJV.l.lmm 


f 1. 1mm 


ln(L) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


J123616.08-f621514.1 


11.10+^30 


901 


4.84+?-5^ 


0.00 


0.001 


-27.94 


J123635.86-I-620707.8 


0.23+n?? 

— U. 1 i 


12 


0.79+o'^? 


0.91 


0.312 


-11.40 


J123711. 32-1-621331.1 




468 


0.72+0:13 


0.01 


0.005 


-40.04 


J123711. 98-1-621325.8 


3.64+'j'-65 


305 




0.16 


0.059 


-9.59 


J123716.63-I-621733.4 


13.00+^5^90 


8825 


17.30+0f3 


0.23 


0.131 


-2868.31 


J033158. 25-274458. 8 


0.67+0-J2 


32 


0.07+003 


0.05 


0.017 


-2226.84 


J033204.48-274643.3 




7 


1.05+0:35 


0.44 


0.021 


-24.87 


J033207.12-275128.6 


0.15in« 


6 


0.04+0:01 


0.25 


0.010 


-23.65 


J033209. 71-274249.0 


, qn + 0.02 


57 


i.87+o:oi 


0.10 


0.070 


-135568.62 


J033212. 23-274620. 9 


1.43+0-^6 


59 


2.94+0:02 


0.11 


0.038 


-3350.01 


J033215. 32-275037.6 


2.04+0:34 


96 


1.10+0:09 


0.01 


0.004 


-1280.28 


J033222. 17-274811. 6 


3.39lo:^? 


269 


3.08+0:^1 


0.18 


0.098 


-60.08 


J033222. 56-274815.0 


r, qn + 0.51 


105 


1.48+0:^6 


0.04 


0.022 


-62.72 


J033235. 18-275215. 7 




3 


0.18+0:01 


0.02 


0.000 


-302.13 


J033238. 01-274401. 2 




167 


0.13+0:02 


0.10 


0.024 


-111.65 


J033244.02-274635.9 


r, c:q+0.92 
z.oo_Q 35 


144 


4 7p;+0.44 


0.00 


0.000 


-89.91 


J033246. 83-275120. 9 




24 


0.30lHi 


0.00 


0.000 


-4001.74 


J033302. 94-275146. 9 


10.70t^;|0 


741 


17.39lJ:g 


0.00 


0.001 


-61.65 


J095905. 05-1-022156.4 


il.4U_3 3Q 


942 


4.08_Q gg 


0.01 


0.009 


-4.53 


J095929. 70-1-021706.4 


9 07+0.94 


133 


Q7+0-21 


0.06 


0.027 


-2.33 


J095945. 15-1-023021.1 


1 11+0.23 
^•J^^-0.01 


49 


1 -lQ+0.04 


0.04 


0.007 


-588.57 


J095953. 85-1-021853.6 




46 


n 70+0.15 


0.00 


0.000 


-222.44 


J095959. 96-1-020633.1 


0.821°:?^ 


65 


36+0 05 


0.51 


0.085 


-174.35 


J100003. 73-1-020206.4 


0.7610-^6 


29 


2 65+0-58 


0.75 


0.212 


-18.55 


J100006. 11-1-015239.2 


l-»U_g og 


78 


9 91+0.57 

^•^^-0.56 


0.21 


0.063 


-11.98 


J100006.55-f023259.3 


1 17+0.21 


54 




0.00 


0.000 


-4319.74 


J100033. 61-1-014902.0 


1 47+0.31 
-0.03 


66 




0.03 


0.005 


-323.60 


J100055. 34-1-023441.1 


Q ^0+0.69 

».00_Q 73 


576 


40.84lJ:98 


0.47 


0.076 


-18184.24 


J100107.46-I-015718.1 


1 40+0-31 


63 


1 C9+O.86 


0.12 


0.028 


-93.88 


J100116.15-I-023606.9 


ap;+0.71 

O.Od_Q g3 


178 


90+1.68 

"•^■^-i.ge 


0.00 


0.000 


-1096.50 


J100139. 73-1-022548.5 


ni+0 02 


9 


04+0-01 


0.32 


0.001 


-187114.26 



limits. Correcting for the under-prediction, we can expect 
the luminosities to be ~2-3 times higher. As a result, the 
FIR-derived SFRs and the associated SFR-derived X-ray lu- 
minosities for J033215.32-275037.6 and J033207.12-275128.6 
are more likely to be in line with their observed X-ray lumi- 
nosities. The remaining starburst-candidate X-ray sources 
(J033158. 25-274458.8 and J100139. 734-022548. 5) are stiU 
ambiguous as their X-ray derived SFRs are ~5-8x higher 
than their IR counterparts; however, the poor fits to these 
sources prevents an accurate measurement of their FIR lumi- 
nosity and SFR, hindering our interpretation. Nevertheless, 
it remains plausible that at least ~6 percent of our X-ray- 
detected SMGs are starburst-dominated in both the IR and 
X-ray with little (if any) emission due to an AGN. 

It is possible to account for the missing 1.1 mm flux 
if we relax the constraints on many of the fit parame- 
ters. For instance, the SB models can provide a better fit 
if we relax the redshift prior. Similarly, if the X-ray lu- 
minosity constraint is removed then the AGN templates 



can account for the remaining 1.1 mm flux with signifi- 
cantly more dust (as set by Ay and R). These fits, however, 
are completely unphysical either due to inaccurate redshifts 
(A2: >0.5) or X-ray luminosity (unconstrained AGN tem- 
plates predict orders of magnitude higher X-ray luminosi- 
ties, see also Figure [Sj. Instead, these additional fits suggest 
that an additional, possibly extended, dust distribution may 
be required. Similar modifications have been suggested for 
other SED templates in order to provide com plete fits to 



other SMGs and millimetre-detected QSOs (e.g. Pope et al 



|200S ; iMarti'nez -Sansigrc et al.' 2009; Rowan -Robinson et al 
[2010 . Unlike (Rowan- Robinson ct al. (201^), however, we 
find that a diffuse 'cirrus' c o mpon ent as described by 
lEfstathiou fc Rowan- RobinsonI (|2003l ) is not sufficient for 
the additional dust distribution and does not improve the 
quality of our fits. 
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Table 6. Best-fit SED parameters using only the SB models. The SED derived X-ray luminosity is left as a free parameter and provides 
no additional constraint to the SED fitting. Errors are given at the Ict confidence level after marginalizing over all other free parameters 
in the fitted templates. Col.(l): Chandra Source ID. Col. (2), (3), and (4): SB template normalization, age, and optical depth. Col. (5): 
Total model-derived, rest-frame bolometric IR luminosity from ~0.001-1500^m. Col. (6): SFR derived from Eqn. 2 (§ 3.2.1) and Cols. 2 
&; 3. Col. (7): SED derived X-ray luminosity. Col. (8): Fractional contribution of model to observed 1.1mm flux. Col. (9): ln(L) of best-flt 
parameters. 
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4 DISCUSSION 

Across a ~1.2 square degree area of the sky, we have ana- 
lyzed the X-ray spectral and NIR-to-radio SED properties of 
45 X-ray-detected AzTEC sources for evidence of AGN and 
starburst activity. Our full sample is limited by the number 
of available redshifts, leaving a subset (32/45) of sources. 
Within GOODS-N and GOODS-S, this subset of AzTEC/X- 
ray sources typically have high levels of dust obscuration 
(Nh > IQ-^-^ cm~^) and are generally associated with AGN 
activity, while their NIR-to-radio SEDs imply that the IR 
and bolometric output are almost completely dominated by 
star formation. Though we do go deeper in the 4Ms GOODS- 

5 field and find fainter potential X-ray counterparts, we do 
not find any evidence for significantly higher amounts of dust 
obscuration compared to the 2Ms GOODS-N and initial 2Ms 
GOODS-S. Considering the relative uncertainties in the Lx- 
SFR relation and under-prediction of the 1.1 mm flux for 
many of our models, a small portion (~6-13 percent) of our 



X-ray-identified SMGs are likely to be completely dominated 
by starburst emission in both the X-ray and NIR-to-radio 
with the remaining majority powered almost exclusively by 
an AGN in the X-ray and starburst in the NIR-to-radio (see 
§ 3.1.3 and § 3.2.2). Here, we explore the implications of our 
X-ray modeling and SED fitting in the context of emission 
at 1.1 mm and previous (sub-)mm/X-ray studies. 

4.1 Origin of 1.1 mm Emission 

As stated in § 3.2, (sub-)millimetre emission from our 
AzTEC/X-ray sources results from dust heated to T~30K. 
Based on our SED fitting (§ 3.2.2), the observed NIR- 
to-mm luminosity is generally dominated by the star- 
burst with little contribution from an AGN (see Figure |3]) . 
These fits predict dust temperatures on the order of ~30- 
40K yet generally under-predict the observed 1.1 mm flux 
by >50 percent, suggesting that a cooler, extended dust 
component is present. Further evidence for an additional 
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Figure 3. Observed frame best-fit AGN+SB and SB-only SEDs to AzTEC/X-ray sources. The plots show the template models that lie 
closest to the best-fit parameters determined from the MCMC SED fitting. The AGN and SB models are given by the dotted red and 
dashed blue lines, respectively, with their linear combination shown by the solid black line. Also shown for each source is the redshift 
used in the SED fitting, favoring the spectroscopic redshift where available, and the resulting best-fit ln(L) values. For reference, we have 
included the probability of random association P (Table [TJ for each source. 



dust co mponent has been seen in previous SM G stud- 



(e.g.. IChapman et all l2005l : iPope et al.l 120061 ) and by 



iRowan-Robinson et al. ( 2010l ) in Herschel SPIRE sources 



although they sugge st that it can be accounted for w i th the 
cirrus templates of lEfstathiou Sz Rowan- RobinsonI (|2003l ') 
which we are unable verify. Fitting of the IR dust peak, 
for which Herschel data are optimized, will provide more 
accurate estimates of the dust temperature and will aid in 
reducing parameter uncertainty, improving the bolometric 
luminosity estimates and providing further insight into the 
nature of the missing dust. 

Given the evidence so far for an additional dust com- 
ponent, one must wonder where the dust resides. The dust 
could simply reside in an extended disk if the starbursting re- 
gion remains localized to the central ~1 kpc. Alternatively, 
the dust may reside in the halo of the SMG, pushed out 
through radiation- or momentum-driven outflows resulting 
from the starburst region(s) and/or the central AGN (e.g.. 



lOppenheimer fc Dave|[200^ : [z^ et al. I I2OIII ). Atypical GMC 
in a z = 2 starburst galajcy can reach vel ocities of ~300 km 
s~^ iMurrav. Menard fc Thompsonll201ll '). which will spread 
its gas and dust as far as ~15 kpc from the galaxy cen- 
ter during a ~50 Myr starburst active phase. Similar out- 
flows reaching ~1000 km s~^ have been observed in local 
ULIRGS and have been shown to account for as much as 
20 percent of the total molecular gas mass, on the order of 
10^ M0, which are ea sily produced through st arbursts with 
SFR> 100 Mq yr-i (|Cliung et al.ll2009l . I2OIII ). The spatial 
scales predicted for these outflow regions are consistent with 
the radii predicted by the AGN templates and high resolu- 
tion imaging of SMGs using the TRAM Plateau de Bure in- 
terferometer (PdBI) and Submillimetre Array (SMA) (e.g., 
iTacconi et al]|2006l . [20081 : lYounger et al.ll2008l ). which show 
typical size scales of ~2-8 kpc. The molecular gas will not 
survive long due to lack of self-shielding, which would allow 
the dust to inhabit a larger volume than that traced by tra- 
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Figure 3 - continued 



ditional molecular gas measurements. Though the majority 
of the mass will still be contained within the central region, 
the extended dust will quickly cool to the background tem- 
perature and will likely produce a temperature gradient as 
distanc e from the central reg ion increases (see, for example, 
fig. 5 of lYounger et al]|2008l ), which may contribute signifi- 
cantly to the (sub-)mm emission. This scenario agrees with 
the recent EVLA observations of Jvison et al. (2 010. 2oTl|) 
SMGs and of background quasars ( Menard et all [201 ol ). 



A possible alternative to the extended cold dust model 
is that the missing 1.1 mm fiux results from either false de- 
tections or source blending. We may readily remove false de- 
tections as the false association rate for the X-ray-identified 
AzTEC sources is r^^S-G percent (§ 2.3.1) whereas the ma- 
jority of sources under-predict the 1.1 mm flux. Similarly, 
previous sub-mm studies suggest blending can o ccur in r^20- 
25 p ercent of sub-mm detected sources (see IScott et al.l 
|2012| . and references therein) so that while blending is likely 
to occur, it is unlikely to the primary cause for the flux 
discrepency. Unfortunately, it is not possible to de-blend 
sources using current Spitzer MIPS, Herschel PACS or 



VLA radio data without a priori knowledge of the intrin- 
sic sources. Only through high resolution imaging and kine- 
matics with ALMA, LMT and future (sub-)mm telescopes 
may we be able to de- blend potential offenders and/or make 
direct confirmation of an outflow-produced, extended cold 
dust distribution. 

The question still remaining is how the AGN emission, 
as indicated by the high X-ray detection rate and the X- 
ray spectra, relates to the sub-mm observations. While AGN 
models are favored in the X-ray spectral fitting, the sub-mm 
emission is, in fact, unlikely to result solely from an AGN. 
As shown in Figure [5l the X-ray priors prevent any sig- 
nificant contribution from the AGN templates. Even when 
relaxed, the AGN models still show poor fits to the mid- 
IR, relative to their observed fiuxes and uncertainties, and 
sub-mm data, never mind the unphysical X-ray luminos- 
ity predicted. In a m erger-driven formation scenario (e.g. 
iNaravanan et al.ll2010l '). gas from the colliding systems gives 
rise to an increase in star formation, resulting in a sub-mm- 
bright phase. Shortly after the sub-mm-bright phase and 
final coalescence, the central black hole may undergo the 
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Figure 4. Smoothed, marginalized likelihood distribution of accepted parameter steps for J123616. 08+621514.1 (top) and J033215.32- 
275037.6 (bottom), using the AGN+SB templates. A description of each parameter in the AGN and SB templates is given in § 3.2.1. The 
location of the maximum likelihood value has been marked by the large 'X'. Contours are drawn at the 68 percent (solid) and 90 percent 
(dashed) confidence levels. The likelihood distributions show that while there are no large apparent correlations between parameters, the 
constraints on some parameters are rather poor (particularly for AGN R and SB Tu). AGN L and SB normalization are the most well 
constrained due to the inclusion of the X-ray luminosity prior. 
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Figure 5. AGN only fits for J03222.56-274815.0 using tfie X-ray 
priors (top) and without (bottom). As seen in Figure [S] the AGN 
component is unable to contribute any more to the NIR-to-radio 
SED as its vertical scaling, i.e. luminosity, is set by the X-ray 
priors. Without these constraints, the models are able to account 
for some of the mid-IR emission, still missing the bulk of the sub- 
mm flux, but would predict tremendous X-ray luminosities; here 
the best fit AGN L is ~ 10^^'^^ Lq which translates to an X-ray 
luminosity of ~ 5.9 X 10^* ergs s~^, many times over what is 
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Figure 6. Histogram of effective power-law indexes (T^ff) for 
the AzTEC/X-ray sources, given by the filled histogram, for all 
sources (top) and only those with radio counterparts (bottom). 
For comparison, we also include the FEff distributions from A05b 
(back-hashed region) and LNPSIO (forward-hashed region). 



actually observed. 



bulk of its growth, producing an AGN which may then aid 
in shutting off the star formation through feedback, leaving 
the final system as a quasar or dusty AGN-powered ULIRG. 
Given the high X-ray column densities we derived for our 
AzTEC/X-ray sample, it is likely that these sources repre- 
sent the early growth phase of the AGN. Combined with 
the starburst-dominated NIR-to-radio SEDs and expected 
short timescale of the sub- mm- bright phase (< 50 Myr), the 
X-ray-identified sources may be SMGs caught in their tran- 
sitionary period between peak star formation and peak black 
hole growth. This transition scenario is consistent with the 
fact that the average 1.1 mm fluxes and 2-10 keV count rates 
of the X-ray-identified sources are both below the average of 
the overall SMG and X-ray sample. The remaining X-ray un- 
detected SMGs could result from a starburst triggered dur- 
ing the first passing of merging syst ems or rapid, short-live d 
mergers similar to those found by iGhapman etal] (|2009h . 
However, we can not rule out the possibility that the X-ray- 
dim SMGs could result fr om a moderately continuous gas 
in-fall (see, for example, iDave et al.l |201(]| ) or very young 



starbursts with Compton-thick AGNs (e.g., A05a,b). The X- 
ray-detected SMGs are unlikely to be produced by such con- 
tinuous in-fall given the starburst timescales from our SED 
modeling; accretion-driven models pred ict that the sub-mm 
bright phase may last for ~0.1-1 Gyr (jPardal et all 120011 ). 
One other possibility given the expected high SFRs for 
SMGs is that the central AGN are likely time- variable (see 
[Alexander fc Hickoxll201^ . and references therein) . AGN can 
switch between being 'on' or 'off' on timescales of <1 yr and 
cause large variations in their observed luminosities and ab- 
sorbing column densities, which will affect the probability 
of detecting an AGN associated with an SMG. It is un- 
known how this AGN time-variability scenario will influ- 
ence the SED of SMGs though we expect any contribution 
to be small given the already low AGN contribution rate. 
Further evolutionary simulations and observations aimed at 
spatially resolving SMGs will provide the tools necessary to 
classify SMGs under the appropriate formation and evolu- 
tionary scenario. 
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Figure 7. Reproduction of fig. 2b from A05b including the A05b, 
LNPSIO and AzTEC/X-ray samples. Tlie fluxes of A05b have 
been converted to 2.0-10.0 keV l uminosities assuming a photon 
index of F = 1.8 and eqn. 1 of [Alexander et al. 1 ||2003|') . Radio 
luminosi ties are calculated from the radio fluxes in Table [2l and 
eqn 2. of lAlexander et al.l ||2003| '). Also plotted is the P04 SFR-X- 
ra y lumino s ity re lation (solid line), using the SFR-radio relation 
of lCondonl ll 19921) to convert SFR to 1.4 GHz luminosity, with a 
20 percent statistical error given by the dotted lines. Some of the 
A05b starburst sources only have 3cr upper limits for their X-ray 
luminosities and are shown with arrows indicating such. 



Figure 8. Reproduction of fig. 8 from A05b including the A05b, 
LNPSIO and AzTEC/X-ray samples. The X-ray fluxes of A05b 
are converted to 2.0-10.0 keV using F = 1.8. FIR luminosities 
are derived from the radio lu minosities of Fig. [7] using a radio 
to FI R correlation of q=2.35 l|Helou. Soifer &: Rowan-Robinsoa 
Il985h . The over-plotted lines represent ratios of constant X-ray 
versus FIR luminosity for the A05b starburst = 10~*) and 

AGN (|-^^ =0.004) sources, and the average luminosity ratio for 

quasars studied bv lElvis et"al] {l99i) {-^^=0.05). 



4.2 Comparison with Previous Studies 

For our AzTEC/X-ray sample, the AGN detection rate is 
~14 percent between all three fields. However, the shal- 
low X-ray depth of COSMOS, potentially compounded by 
our more stringent detection criteria, prevents confirma- 
tion of the most heavily obscured AGNs whic h may con- 
tribute significantly to t he sub- mm emission (|Lutz et al.l 
I2OIOI : iHill fc Shanks! bOllh . Excluding COSMOS, the AGN 
detection rate increases to ^ 28 percent, consistent with 
previous X-ray/SMG studies (> SSIjq percent, 29±8 per- 
cent and < 26 ± 9 for A05a,b, LNPSIO and GRCll, respec- 
tively) while avoiding potential biases due to prior coun- 
terpart identification and achieving better source statistics 
via larger sky coverage. Similar to LNPSIO, we also find 
evidence that ~6-13 percent of our X-ray sources are po- 
tentially HMXBs associated with high star formation rates. 
However, many of the starburst powered SCUBA-detected 
sources of LNPSIO, and by extension A05b, are missing from 
our sample. While the our X-ray data for GOODS-N is es- 
sentially the same as that used in A05b and LNPSIO, it is not 
surprising for diff erences to exi s t betw een the AzTEC and 
SCUBA catalogs. IChapin et all (|2009| ) suggests that such a 
discrepancy results from instrument and measurement cali- 
bration uncertainty as well as intrinsic spread in host prop- 
erties (namely dust temperature and emissivity). In fact, 
for a SCUBA source to be detected by AzTEC at > 3.5a in 
GOODS-N (where the AzTEC rms is ~1.3 mjy/beam, see 
§ 2.1), its effective 850/im flux would need to be > 8.19 
mjy, higher than the typical 850/im flux for sources in 
LN PSIO. This estimat e assumes an R=S85o/Si.i value of 
1.8 llChapin et al.ll2()09l) and that 'flux boosting' (see, for ex- 
ample, Austermann et aLll2010l : IScott et al.ll2O10l ) effects the 



850/im and 1.1mm observations equally. Completeness of the 
(sub-)mm observations may also contribute to this discrep- 
ancy; at ^4mJy, th e AzTEC map is ~60 percent complete 
(|Perera et al.ll200^ . Of course, there is always the issue of 
false identifications and mismatching of sources as well as 
prior counterpart bias (see LNPSIO) which, while the ex- 
pected number of such occurrences are small (see § 2), may 
still lead to a decrease in the number of starburst-dominated 
X-ray sources in our sample. 

In Figure [6] we show the range of effective photon in- 
dexes Feh for our AzTEC/X-ray sample (see § 3.1.1, Ta- 
ble [3]) in relation to the samples of A05b and LNPSIO. Us- 
ing the Mann- Whitney (MW) U-test, we find that the prob- 
ability that our AzTEC/X-ray sources are consistent with 
being drawn from the samples of A05b and LNPSIO are 0.02 
and 0.14, respectively. If we limit our sample to AzTEC/X- 
ray sources with radio detections then the MW probabilities 
become 0.07 and 0.17 for the A05b and LNPSIO samples, 
respectively. Since the errors on FEff are known, we further 
estimate the intrinsic mean and variance of the samples by 
constructing 1000 Monte Carlo realizations of the Feh dis- 
tributions. The resulting intrinsic mean value of Feh for the 
AzTEC/X-ray, A05b and LNPSIO samples are 1.14±0.09 
(Icr), 0.60±0.10 and 1.44±0.16, respectively; including only 
the radio-detected AzTEC/X-ray sources results in an in- 
trinsic mean of 1.05±0.08. These results imply a strong sta- 
tistical difference between the AzTEC/X-ray and A05b sam- 
ples (at >3cr), while the AzTEC/X-ray and LNPSIO samples 
have consistent means values of FEflf. 

Despite the differences in FEff, the methods of analy- 
sis in A05b produce results consistent with our study. For 
further comparison, we reproduce figures 2b and 8 of A05b, 
which show the L2.0_10.okcV versus Li.4ghz (Figure [7]) and 
L2.0-10.0keV versus Lpir (Figure [8]) relations for the A05b 
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and LNPSIO starburst and AGN systems, including our 
AzTEC/X-ray sources that have radio counterparts. In re- 
producing the A05b figures, we have converted the A05b 0.5- 
8.0 keV fiuxes to 2.0-10.0 keV luminosities assuming a pho - 



ton index of F = 1.8 and eqn. 1 of [Alexander et al.r(|2003l ). 
Radio and FIR luminosities have been determined for our 
sample following the same procedures as A05b to ensure 
compatibility. We caution, however, that the radio-FIR cor- 
relation u sed to derive the FIR luminosities from the radio 
emission l|Helou. Soifer fc Rowan-RobinsonI Il985l ) assumes 
emission purely from star for mation and could be misle ading 
if the AGN is radio-loud (e.g. lPonlev et al.ll2005l . l2010h . Fig- 
ures[7]and[8]show that the X-ray emission for the sub-sample 
of radio-identified AzTEC/X-ray sources is higher than one 
would predict from their radio and/or FIR luminosities if 
they resulted purely from star formation, indicating AGN 
activity. However, the FIR luminosities are generally higher 
than expected for typical quasars which suggests significant 
contribution from star formation, again consistent with the 
results from § 3.2 . Alternatively, sources could lie above the 
Hi vis et all (j 19941 1 quasar relation if they are reflection dom- 
ina ted or Compton -thick (e.g. FSC 10214-1-4727 A05b, Arp 
220 llwasawa et ail '2005). This is not likely to affect our anal- 
ysis based on the results from our X-ray spectral modeling 
(§ 3.1); nevertheless, we can not rule out the possibility that 
the faintest X-ray sources may be harboring highly lumi- 
nous, Compton thick AGNs, particularly f or the non-X-ray- 
detected SMG s (e.g . Ilwasawa etHI l2005l : Eutz et al.ll2010l : 
ImU fc Shank3l201lh . 

4.3 Cross-Correlation of AzTEC/X-ray source 
populations 

In addition to examining X-ray-detected SMGs on a source- 
by-source basis, a simple cross-correlation analysis of the 
X-ray and AzTEC source populations ca n identify evolu - 
tionary patter n s bet ween the populations. lAlmainil | 20021) : 
lAlmaini et al.l (|2003l ) were the first to measure the angu- 
lar cross-correlation function (XCF) between SMGs and X- 
ray sources and found significant correlation at large scales, 
leading to the conclusion that the populations both reside 
in similarly massive dark matter ha l os and trace the same 
large scale structure. iHill fc Shanks! (|201ll ) later estimated 
the XCF for LABOCA sources in the ECDFS and while 
found similar evidence for small scale clustering, i.e. re- 
siding in same dark matter halos, fo und no evidence fo r 
large scale clustering, con s istent with iBorvs et al.1 (|2004h . 
iRoche. Dunlop fc Almaiiiil (|2003D measured the XCF of ex- 
tremely red objects (EROs) - which may be the signature of 
massive galaxies that have entered their passive post- AGN 
phase in galaxy evolution - and X-ray sources in the CDFS, 
and again find evidence for significant correlation. Together, 
these results may suggest an evolutionary sequence between 
these three populations, where starburst dominated SMGs 
go through an AGN-bright phase before evolving into pas- 
sive ellipticals or EROs. 

To determine if there is any correlation between the 
AzTEC and Chandra source populations, we apply the two- 
point angular XCF, wax{0), defined as the excess proba- 
bility of finding both an AzTEC source in a solid angle 
5^ A and an X-ray source in a solid angle 5Q.x, with an 
angular separation 6 from each other. This excess proba- 
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Figure 9. The XCF between AzTEC and Chandra source popu- 
lations. Plotted in each panel is the observed XCF and the XCF 
from randomly generated source populations along with the re- 
spective beam-size and search radius for each field. Below our 
adopted search radii, the XCF shows significant signal due to de- 
tected counterparts. The lack of consistent positive correlation in 
COSMOS results from the shallow X-ray depth and correspond- 
ing low source density. 



bility (relative to an uncorrelated distribution) is given by 
5P — pAPx[i + 'WAx{d)]SilA5iix , where pA and px are 
the surface densit ies of AzTEC and X-ray galaxies on the 
sky (|Peebleslll98Gl l. In practice this can be measured from 
galaxy maps by counting the number of SMG/X-ray source 
pairs, binned by their angular separation, and comparing 
to pair counts from random positions. Here, we use the 
cross-correlation adapt ation to the Landy-Szalay estimator 
(|Landv fc Szalavlll993h . which is given by 



WAx{0) = 



DADx{e) - DARx{e) - DxRa{9) + RaRx{9) 



RARx{e) 



(3) 
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Figure 10. The XCF between the AzTEC and Chandra source 
populations (shown in Figure at scales larger than the beam- 
size. The top panel shows the full sample while the bottom shows 
the XCF in the redshift range 1< z <3. Over the three fields, 
there is no significant correlation between the source populations, 
particularly over the typical redshift range of SMGs. Due to sen- 
sitivity variations in the Chandra data, the XCF should not be 
heavily weighed at angular separations of >200. 



where DaDx is the number of SMG/X-ray source pairs, 
DaRx and DxRa are the number of pairs found between 
each galaxy catalog and randomly generated positions of 
sources within each angular separation bin. RaRx is the 
number of pairs found between random positions for each 
galaxy population, generated from the selection function 
and sensitivity distribution of each map. To generate ran- 
dom source distributions for the A zTEC maps, we follow the 
methods o f lWiUiams et al.l (j201lh . For the Chandra random 
catalogs, the exposure maps are relatively uniform (ignoring 
effects due to CCD gaps and edge overlapping in COSMOS 
as they are generally small) such that we may produce the 
random catalogs by simply randomly generating positions 
within the overlapping coverage region of the Chandra and 
AzTEC maps. Note, however, that this does not take into 
account the sensitivity variations (mostly due to PSF degra- 
dation) as a function of off-axis distance; the XCF may thus 
be incorrect at scales larger than ~200". The overlapping 
observations in COSMOS helps to smooth the telescope re- 
sponse, allowing for a more accurate XCF at larger scales. 
The resulting XCF for each field, as well as the ex- 



pected XCF from completely random distributions, is shown 
in Figure [9l where the errors are estimated from a Poisso- 
nian distribution given the number of AzTEC/X-ray pairs 
in each angular bin. The expectation from random distribu- 
tions is estimated by averaging the XCF of 100 AzTEC and 
X-ray random distributions described above, which have the 
same properties (area and source density) as the observed 
maps. In the case of the random expectation, the errors cor- 
respond to the standard deviation of the XCF from each of 
the individual random distributions. At small scales, there 
is significant positive correlation in the observed XCF due 
to identified counterparts (see also lHill fc Shank^l201ll ): this 
effect is diluted in COSMOS due to its shallow X-ray depth 
and thus low source density compared to either GOODS 
field. However, since the AzTEC source positions are not 
well known on scales smaller than the beam-size, we choose 
to limit our XCF analysis to the large scale clustering. Fig- 
ure [TD] shows the same XCF combined with their weighted 
average for scales larger than 28", the beam-size of AzTEC 
on ASTE, though we caution against heavy interpretation 
at scales larger than ~200" as previously mentioned. 

Across the three fields, we find no evidence for any large 
scale correlation signal; any apparent correlation or anti- 
correlation seen in ind ividual fields is det ecte d at < la confi- 
d ence , consistent with lBorvs et al.] l|2004 ) and lHill fc Shank^ 
(2011). The large area covered by our sample (~1.2 square 
degrees) aids in mitigating the effects of cosmic variance, 
which is the likely cause of variation betwee n fields and ma y 
affect the positive correlation signal found in lAlmainill|2002h . 
It is possible that the lack of any correlation signal in our 
data may be the result of dilution given the wide and differ- 
ing redshift distributions of the X-ray and sub-mm sources. 
In an attempt to improve the cross-correlation signal, we 
have run the same analysis by limiting the X-ray sources to 
the redshift range of 1 < z < 3 where the X-ray redshift 
distribution shows significant overlap with the sub-mm dis- 
tribution. If there is any cross-correlation between the two 
samples, it should be maximized here. Due to the small num- 
ber of X-ray sources with available redshifts in GOODS-N 
(49 out of the original 397), we excluded this field from the 
XCF in the 1< z <3 redshift range. The XCF using this 
redshift-limited subset for GOODS-S and COSMOS is sta- 
tistically identical to the result we measured using the entire 
set of X-ray sources, i.e. no evidence for a correlation. 

The lack of a significant correlation between the X-ray 
and AzTEC source populations at large scales may suggest 
that SMGs and AGN are not universally related in terms 
of dark matter halo mass and large scale structure. How- 
ever, considering the significant fraction of AzTEC SMGs 
that do have plausible X-ray detections here, it is likely that 
the SMG phenomenon is not governed by a single forma- 
tion and evolution process; rather, the SMG population is a 
"mixed bag" of systems - some undergoing major mergers 
concurrent with the b uild-up of massive black holes (e.g., 
[Narayanan et all |2010| ) and others signaling a short-lived 
phas e of intense star-form ation in more normal galaxies 
(e.g.. lChapman et a l."2009''l or e ven qu iescent mass build-up 
from gas in-fall (e.g., Dave ot al.ll2O10l ') (see also § 4.1). Such 
cases are likely tied to the host's intrinsic properties which 
could naturally explain the enhanced sub-mm emission from 
bright, obscured AG Ns as found by iLutz et al.l (|2010l ) and 
iHill fc Shan"k3(|2011^ . However, we caution that limitations 
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in measuring the correlation between these populations can 
also give a null result. For example, the large volume sam- 
pled coupled with the lack of redshift information for the full 
X-ray and SMG catalogs will necessarily dilute the projected 
correlation strength between the two populations, even if 
there is some spatial correlation. The shallow X-ray depth 
of COSMOS will further dilute any correlation signal by pri- 
marily detecting bright AGN that are likely well past their 
starburst phase. Observations of SMGs in the near-future 
with ALMA and the LMT geared towards measuring their 
redshifts and obtaining high-resolution imaging of their dust 
and gas will greatly aid in the development and fine tuning 
of formation and evolution scenarios for this population. 



5 SUMMARY 

We have presented a detailed analysis of the X-ray prop- 
erties of AzTEC 1.1 mm sources found in the GOODS-N, 
GOODS-S and COSMOS fields. Thanks to deep (~2-4 Ms) 
Chandra observations, we find X-ray counterparts to ~14 
percent of the 1.1 mm sources across all three fields, in- 
creasing to ^^28 percent if we exclude COSMOS due to its 
shallower X-ray data. From our modeling of the X-ray spec- 
tra and NIR-to-radio SEDs, we conclude that AzTEC/X-ray 
sources arc all starburst-dominated in the IR, with SFRs on 
the order of 100—1000 MQyr^^, whereas an AGN component 
is needed in order to explain the observed X-ray luminosi- 
ties for the majority of our sources. In ~6-13 percent of our 
sample, we find evidence for X-ray emission consistent with 
high SFRs, after accounting for the relative uncertainties 
in the Lx-SFR relations and the typical under-prediction 
of the 1.1 mm flux in our SED modeling. The AGNs typi- 
cally appear obscured in the X-ray band, with neutral hy- 
drogen column densities in excess of 10^^ cm~^. These re- 
sults are consistent with other SMG/X-ray studies. Overall, 
the AGN templates contribute very little (<10 percent) to 
both the bolometric luminosity and 1.1mm flux. At 1.1 mm 
in particular, the AGN-I-SB models typically under-predict 
the observed fluxes, which indicates that either a cooler, ex- 
tended dust component is required to fully recover the NIR- 
to-radio SED or that the sources are blended. We suggest 
that this missing dust could result from radiation- and/or 
momentum-driven outflows caused by the starburst /AGN 
regions, which pushes the dust out into the halo where it 
cools rapidly and, although it accounts for a small fraction 
of the total dust mass, may contribute significantly to the 
1.1 mm emission. 

The high AGN identification rate in these AzTEC 
SMGs is particularly interesting in regards to SMG forma- 
tion and evolution scenarios. Following a merger-driven sce- 
nario, the X-ray identified sources could represent the transi- 
tional period between starburst and AGN dominant phases. 
However, the lack of a significant correlation at large scales 
between all X-ray sources and SMGs in these fields suggests 
that not all SMGs will evolve to possess an AGN and, simi- 
larly, that not all AGN evolve from a sub-mm bright phase. 
This suggests heterogeneity in the formation/evolution of 
SMGs, possibly due to either intrinsic source properties, i.e. 
amount of obscuration, or even multiple formation scenar- 
ios. With future analyses aimed at source evolution as a 
function of redshift, combined with a more comprehensive 



redshift catalog for SMGs (one of the goals for the upcoming 
LMT), we will be able to determine the AGN fraction and 
contribution to greater certainty, allowing for investigating 
how SMGs form and evolve into the galaxies we see in the 
local Universe. 
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